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We investigate the consequences of introducing a set of exotic doublet leptons which couple 
to the standard model leptons in a minimal way. Through these additional gauge invariant and 
renormalizable coupling terms, new sources of tree-level flavor changing currents are induced via 
mixing. In this work, we derive constraints on the parameters that govern the couplings to the 
exotic doublets by invoking the current low-energy experimental data on processes such as leptonic 
Z decays, £ — ► 3£' , £ — >• i'^y, and \i-e conversion in atomic nuclei. Moreover, we have analyzed the role 
these doublets play on the lepton anomalous magnetic moments, and found that their contribution 
" is negligible. 
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PLh ' I. INTRODUCTION 
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It is well-known that the phenomenon of neutrino oscillations |l| have motivated an extensive study on models with 
nonzero neutrino masses and lepton flavor violation (LFV). Central to all these investigations is the introduction of 
' (— i ' 1 new interactions (and most likely new particles too) to the minimal standard model (SM). While there can be many 
O ' y different ways of inventing and restricting the new physics, recently in 0, a fairly economical approach based on only 
' ■ SM gauge invariance, renormalizability and the concept of "minimal coupling" was considered for the lepton sector. 

In that work, a generic minimal coupling between SM particles and some exotic field was defined to have the form 



Oh 



O 



^exotic (SM particle) • (SM particle) • (exotic particle) , (1) 



CNI , where Exotic denotes the coupling strength and the "exotic particle" can either be a scalar boson, a fermion or a 
vector boson. All particles in (JTJ) are assumed to be uncolored in the SU(3) C sense because we only wish to extend 
the lepton sector. These minimal interactions are of interest because they are relatively simple and their collider 
• signatures may be detected at the LHC in the near future under favorable conditions Q. 

Schematically, there are five distinct types of interaction with the SM fields allowed by ((TJ) and the requirement of 
renormalizability 1 , 



(i) Ll x Ll x [new] , (ii) Ll x £ R x [new] , (hi) £r,x IrX [new] , 

(iv) Ll x (jy x [new] , (v) Ir x <j) x [new] . (2) 



where Ll = {vl,^l) T is the left-handed (LH) lepton doublet, £r is the right-handed (RH) lepton singlet, and 
4> = ((f> + , 4>°) T denotes the SM Higgs doublet. The imposition of Lorentz and SM gauge symmetries will then result 
H ■ in only 13 types of exotic multiplets that are minimally coupled to the SM particles (see Table 

As pointed out in 0, it seems that all but two types of the new particles induced by this setup has already been 
heavily analyzed because of various motivations. The two that are rarely discussed arc the exotic lepton triplets 
Eb,,l = (El,r>E l ~ r,El~r) T an d the doublets Ll,r = {L^ R , L~[~r) t with SM quantum numbers (1,3,-1) and 
(1, 2, —3/2) respectively. Whereas the study of the triplets Erl formed the backbone of @, we shall concentrate on 
the last remaining possibility — doublets Ll.r, in this work. 

The various implications of introducing Ll.r. to the SM are studied in the subsequent sections. Processes such as Z 
decays (Sec. IIIip. LFV decays: £ — > 3£' (Sec. ITV)) . £ — > i'j (Sec. IV]) . and /i-e conversion in atomic nuclei (Sec. IVI[) are 
considered with the aim to derive constraints on the relevant new physics parameters using low-energy experimental 
data (Sec. IVTip . The contributions from Ll,r to the lepton anomalous magnetic moments will be investigated in 
Sec. I Villi Finally, in Sec. IIX1 we will make some brief comments regarding their collider phenomenologies. 
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1 We do not consider terms such as (SM Higgs) (SM Higgs) (new boson) where no leptons of any type is present. 
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TABLE I: Summary of the 13 types of exotic multiplets induced by the 5 general types of minimal couplings: (i) LlxLlX [new], 
(ii) Ll x £rX [new], (iii) £rX £rX [new], (iv) Ll x 4>x [new], and (v) Ir X cj>X [new]. Hypercharges are defined with Q — I3+Y. 



II. MODEL WITH EXOTIC LEPTON DOUBLETS, L l ,r 

We begin by writing down the framework and explaining our notations for our minimally extended SM with exotic 
lepton doublets. In this model, we have two new sets (LH and RH) of lepton doublets Ll,r, all having hypercharge 
—3/2 (where Q = I3 + Y). In matrix form, they are given by 



Li 



L\ 



L, 



(1,2,-3/2) 



where Lt, and L R are independent fields. 2 The interaction Lagrangian of interest is 3 

I = L L i$L L + L R 0L R - [Z L ML R + L l Y^ c £r + L L Y^t R + L L Y^ C v R 
where <jf = (0°*, — 4>~) T , and the covariant derivative is defined as, 



V2 




cos ( 



■t(h - sin 2 9 W Q) + ie 



(e>0) 



(3) 



(4) 



(5) 



with Q and J3 being the operators for electric charge and the 3rd component of isospin respectively. In the 
Yukawa term involving Y< defines the minimal coupling between £ R and Ll while M sets the energy scale of the 
new physics. Notice that one cannot have a similar type of minimal coupling between L R and any other SM leptons 

because of SM gauge invariance. But its effects on the SM sector can enter indirectly via the mass term LlML r . 
Writing out all the relevant interactions in (|4]), we have 



c + C SM = c 



<w 



C z + C V 



H 



(0) 



2 Because of the identical transformation properties for LH and RH fields, chiral anomalies cancel automatically in this setup. Furthermore, 
by introducing a pair of these, we have maintained an even number of doublets in the overall model, and hence, avoiding any issues 
with global SU(2) anomalies 1231 . 

3 We have also introduced three RH neutrino fields, ur, so that neutrinos can have a Dirac mass. This is done so because we know that 
neutrinos are massive. However, to avoid making the model too complicated and the risk of masking the effects from the exotic doublets, 
we have opted not to include a Majorana mass term for simplicity. But we shall briefly comment on the effects of the Majorana mass 
term at the end of this section. For a full discussion though the readers may refer to the work of lilt |2:| . 
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where 

c w = 
c z = 

^mass 

c H = 



9 

9 



cos 9„ 



1 



V L %V L + 



sin 2 6 W I L R $L 1 



+ h.c. , 

f sin 2 6 W IrZIr 
2 sin 2 9 



2 sin 2 t 



$L R 



-L~ L ML- R -L L ML R 



7=L L Y^ Ir — ih"fn^R. — VLiTiDVR + h.c. 



V2 



7=L T Y^IrH l^m^RH v^ra^VRH + h.c. 

V2 v v 



(7) 

, (8) 
(9) 
(10) 



In getting ([9]) and (fT0|) . we have written (f> = ((f> + ,(j) ) T = {<p + , (v + H + irj)/\/2) T , where v is the Higgs vacuum 
expectation value (chosen to be real), r\ and </) ± are the would-be Goldstone bosons. Also, we have defined mg = 
vYi/^/2 and mu = vY y /^/2. 

The mixing between the SM charged leptons and components of the exotic doublet may be readily derived if we 
write the relevant terms in £ mass in matrix form: 

7TL£ 

(11) 

where we have used the fact that ifji^ = ip&Pl l° r an y fermion field '01,2 ■ One may define the following unitary 
transformations to bring all fields into their mass cigenbasis: 






V L ,R 




lL,R 



J L,R, 



= Ul,r 



"L,R 



L 



L,R, 



^L,R 



= In 



^L,R 

S l ~r7l) c 



(12) 



where the subscript m indicates the mass basis. In ([12")) , Vl,r, Ul,r and /„ have dimensions 6 x 6, (3 + n) x (3 + n) 
and n x n respectively, with n being the number of generations of the exotic Ll,r. fields added to the model. 

Without loss of generality, one can choose to work in the basis where me and M are real and diagonal. As a result, 
/„ is in fact the identity matrix. Furthermore, to make the notations less cluttered, we absorb the neutrino right 
diagonalization matrix XJ v r into mo- In other words, we set mo = m' D U^R, where Ul L m' D U^R = ?ri^ ag and U v l is 



the neutrino left diagonalization matrix. With these conventions, and to O (v 2 M , 
are given by 



the transformation matrices 




U L = 



Vr 



1 vmiYM-i/V^ 
-vM-^me/Vl 1 ; 



Ur = 



1 

.0 u, 



vL. 



(13) 
(14) 



where 



A = — Y M~ 2 Y^ , and 
2 



A' = — M^Y^YM' 1 
2 



(15) 



are 3x3 and nxn matrices in flavor space respectively. Note that U v l may be identified as the usual neutrino mixing 
matrix, J7 PM ns- 
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The relevant terms in the interaction Lagrangian with respect to the mass eigenbasis arc therefore given by 



C w = 



V2 
9 



Lr 



1 L-) ^ + [Pl 9 7l + Pr9Vr] 



\ J i 



Pl9?l 



I L-) Z[PLg$ L + P R gL] 



L— (L—y) tPL9% L 



\ / m 



J L-) H[P L g% L + P R g% R ] 



L~ 



(16) 



(17) 
(18) 



with the new generalized coupling matrices given by (to leading order) 



n w 

9lL 



9lL 



z 

92R 



9lR 



Vl L vUt L meYM- 2 /Vr 



9% 



-vm e YM~ 2 /V2 
1 0, 



9Vn = o 



92R ~ 



-vYM^/V2 
1 - A' , 



1/2 
0, 



92L 



-1/2 + sin 2 9 W mtYM- 2 /\/T 
<M- 2 Yimt/V2 l/2 + sin 2 6» lu , 



9Il = 



-1/2 + 2 sin 2 V 






1/2 -2 sin 2 6 U 



sin 2 9 W + (1/2 -sin 2 9 W )X 



• YM- 1 /2V2 



~uAf" 1 rt/2V2 1/2 + sin 2 6 W - (1 + sin 2 9 W )X' 



Ul L m D /v 







0, 



92R 



m e (2X-l)/v -m e YM^-/VV 
Y\X-l)/V2-M- 2 Y^m 2 e /V2 —v Y^YM -1 /2 j 



92L 



h \t 



(19) 

(20) 
(21) 



Note that there is no need to include gf Rl g^ R and g^ L because they are redundant when the lepton fields are grouped 
in this matrix form. 4 The setup described above will allow us to easily study the new phenomenologies and any 
subsequent constraints of introducing these exotic doublets. 

Any new contributions to tree-level flavor changing currents will be provided by the nonzero off-diagonal entries 
of matrix A. For instance, the presence of A in the upper-left (3 x 3)-block of g^ R is indicative of this fact. But 
despite the introduction of new mixing effects in certain sectors of the theory, one observes that the SM charged 
current interaction remains unaltered at tree-level. This is understandable since the coupling to the W boson is only 
non-trivial for LH particles in the SM while the exotic fields connect to the RH sector exclusively. As a result, this 
also explains the reason A enters only in the P R term of the SM neutral current Lagrangian (in index form) : 



C 



NC 



9 



■Pl-^Vi 



sin 



Pl 



r- ) 5 a p + Pr : - - sin 9 W X a 



1 



(22) 



A related observation regarding the differences between this model and the one studied in [2( is that the modification 
to the value of the Fermi constant as extracted from muon decay experiments (fx — > e + missing energy) will now enter 
at O (A 2 ): 



(Gi 



G|(l+C|A eM | 2 ), 



(23) 



where Gf — V% g 2 /8M^, and G' F denotes the new Fermi constant in the presence of the new physics invoked by the 
exotic doublet L, while C is a numerical factor of O (1). Hence, to a very good approximation, we may simply take 
G' F ~ Gf (ic. A independent) in all our calculations. 



4 We point out in passing that in Q, the corresponding coupling matrices which mix the ordinary charged leptons (£) with the doubly 
charged exotic particles (cf. R in this model) were not investigated. While their omission would not affect the overall conclusions 
reached in 0, they do result in a slight change of the numerical values coming from analyzing the I — > £'7 graphs. It is therefore more 
complete to include them as has been done here for this doublet model (see Sec. IVll. 
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Comments on the neutrino Majorana mass term 

In this subsection, we make a quick digression to comment on the effects one would get if the neutrino Majorana 
mass term, —Mr (vr) c is included in Lagrangian As it is well-known that such Majorana mass term will 

induce mixing between the fields and vr, leading to light and heavy mass eigenstates. The amount of such mixing 
is characterized by the new physics scale Mr (or more precisely the ratio, Y v /Mr). We shall demonstrate below that 
this effect can eventuate in the modification of the interaction Lagrangian that is akin to the role played by A and A' 
from earlier. 

For the following discussion, we shall assume that there are three generations of vr and (without loss of generality) 
that the 3x3 mass matrix Mr is real and diagonal. Upon including the Majorana mass for vr in ((4]), it would be 
more convenient to rewrite the first term in (jlip as 



As usual, one can turn the fields into mass eigenstates via unitary transformations: 
with Vl defined by 

xt i ( 1-e -moMp 1 ] (u vL 0\ ( {1 - e)U vL 





and 



e = V — Y v M R 2 Y} = m D M R 2 rrf D ; e' = V — M R 1 Y} Y v M r x = M^ m D M^ . (27) 

This transformation will lead to the typical type-I seesaw result for light neutrinos: U^Ltn^U^ ~ — m^M^m^. It 

is important to note that e, e' ~ O (M R 2 ) can be comparable to A, A' ~ O (M~ 2 ) if the two new physics scales are 
similar. Hence, this is the reason we have introduced them in definition (|26[) . Physically, the quantity e represents 
the flavor mixing correction to the LH neutrino kinetic energy terms. Such effect is originated from a dim-6 gauge 
invariant operator in the effective Lagrangian as discussed in [231 ] . After spontaneous symmetry breaking, the dim-4 
kinetic terms receives small contribution from this dim-6 operator, which then leads to a non-unitary mixing matrix 
for ordinary leptons. This result is evident from the modified SM interaction Lagrangian for the charged and neutral 
currents after introducing R and v R (with Majorana mass): 

£% c - A=vi^ + p L [([/t i ) fa (^ - e ^)] t p + h.c. , (28) 



C2 ° = c^ {~^ Pl (t ~ ( U lL)i<*e a dU»LhA v 3 +l a 



sin 2 6 W ^ j 5 a p + Pr ( i - sin 2 6 W ) \ a 



(29) 



Because matrix e can be non-diagonal, the quantity inside the square brackets in ([25)1 is non-unitary in general. 
Comparing CJ£ C with ([22)) . it is clear that the model is now more complicated as there are two new effects (e and A) 
entering. In fact, many of the entries in the generalized coupling matrices ([T9")) to (f2"T)) will also get modified with Mr 
or e dependent terms. This, however, should not come as a surprise since Vr is itself an exotic particle, very much 
like the doublet Llm (see Table [J). 

In the light of this, we have opted to omit the additional mixing effects caused by the RH neutrino Majorana terms 
from our analysis. Alternatively, one may think of this as taking the limit Mr 3> M, so that only the A mixing effects 
will be important. 

III. CONSTRAINTS FROM Z DECAYS 

As hinted earlier, the key to any new physics contributions to the electroweak processes can be parametrized 
by the elements of the A matrix for they incorporate all the essential information regarding the exotic doublets L 
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(see (|15p for definition). Therefore, it is useful to study the phcnomenological constraints on the entries of A from 
precision measurements. In this section, we investigate the bounds coming from tree-level Z decays to charged 
leptons: Z — > £ a £/3- The a = (3 cases will place restrictions on A aQ 's whereas for a/ft the off-diagonal entries can 
be constrained. 

Although the limits presented here for the off-diagonal elements will not be as stringent as those obtained from 
other LFV interactions (see Sec. HVl IVl and rVI|) , the constraints for the diagonal elements of A will be useful in the 
later analysis of the anomalous magnetic moments (see Sec. I VIlTj) . 

Calculating the decay rate, T(Z — > £ a £a), using the standard method but with the modified couplings in (|22p . we 
obtain (in the limit of massless final state leptons) 



T{Z -> £ a £ a ) 



G F M\ 



1 



— sin 



Aqc 



(30) 



where W and Mz are the usual Weinberg an gle and Z boson mass respectively. Inserting the decay widths and values 
of the constants obtained from experiments [24[, Eq. (|3T)]) leads to the following constraints for each lepton flavor a: 



'6.1 ±0.3 x 10~ 3 
X aa < \ 7.0 ±0.3 x 10~ 3 
8.0 ±0.3 x 10" 3 



a 
a 
a 



-■ e , 
: M 

: T 



(31) 



Having established the constraints for X aa above, one may further estimate the change to the polarization asymmetry 
of the Z — > 11 decay due to the effects of the exotic doublets. In terms of the new generalized couplings, we have 



\l — (92L)ll (glffjll 

= G&)?i + G&)?i 



1/4 - sin 2 w -2 sin 2 w (l/2- sin 2 w )\ aa + O (A 2 ) 
1/4 - sin 2 w + 2 sin 4 W + 2 sin 2 9 w (l/2 - sin 2 6 w )\ aa + O (A 2 ) ' 



(32) 



where the subscript "11" indicates the appropriate matrix element. Note that the SM prediction for this asymmetry 
can be recovered from (f3"2"j) by setting X aa = 0. At the Z resonance, it is well-known that the forward-backward 
asymmetry for the process e + e~ — > £ + £~ is related to A £ LR simply via 



(33) 



For our exotic doublet model with sin 2 W = 0.231, this quantity is evaluated to about 1.6 x 10~ 2 if X aa = 6.0 x 10 -3 
while it is about 1.7 x 10~ 2 if X aa = 1.0 x 10~ 3 . Comparing this with the experimental best fit (24[, 1.71 ±0. 10 x 10~ 2 , 
one can see that the limits established in (|31[) arc more or less consistent with this precision test. The small deviation 
which appears between the two choices of A QQ demonstrated above indicates that for all calculations, the most 
conservative approach is to interpret (|31j) as 



X aa < 10 3 for all a 



(34) 



Next, we consider the case where a//). The decay rate is given by 

r3 / -, \ 2 



(35) 



Note that in the limit A Q( g — > 0, this rate disappears in accordance with the fact that there is no flavor changing 
neutral currents (FCNC) at tree-level in the SM. Writing this as a branching ratio and keeping only the leading order 
terms, one obtains 



Bx(Z -> £o£b) = F(Z ^ Br(Z -> £~£ a ) , 
V 1 ' T{Z £ a £ a ) V ' 

_ |A a)3 | 2 (4sin 4 ^-4sin 2 W + 1) 
8 sin 4 9,„ - 4 sin 2 m + 1 



Br(Z ->■ l a l a ) 



(36) 
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From this, we can derive the following bounds for |A Qi g|: 5 

|A eM | < 6.6 x 10- 3 , (37) 

|A eT | < 1.6 x 10" 2 , (38) 

|A Mr | < 1.8 x 10~ 2 . (39) 

Because we are working in the basis where M is real and diagonal, A is necessarily hermitian. Therefore, we have 
\\xp\ = |A^ Q |. As a result, all entries of the A matrix can now be constrained. 

But as foreshadowed, we know from past experience that some of the strongest bounds on such new physics would 
come from the lepton flavor violating decays of ordinary charged leptons. So, processes like £ — > 2>£' and £ — > £ ,r y are 
expected to yield even stronger bounds for \X a p\ than those presented in this section. Moreover, we anticipate that 
the strongest limit on |A eA1 | will come from the studies of muon-to-electron (/i-e) conversion in atomic nuclei as it is 
well-known that this process gives rise to a very strong constraint on the fi-e-Z vertex (25|. We shall investigate these 
in following sections. 

IV. TREE-LEVEL I -> W DECAYS 

Assuming three generations of ordinary leptons, there are only three generic types of final lepton states possible for 
a charged lepton decaying into three lighter ones: t pi pip, £ a £p1p and £ a £p£p, where /3 7^ a 7^ a, with a denoting the 
flavor of the decaying lepton. In theory, the mediating particle here can either be the gauge boson Z or the Higgs 
boson H. However, since the amplitude associated with the Higgs is suppressed by a factor of m 2 /M^ (where m a 
and Mh denote the lepton and Higgs masses respectively), it is safe to neglect their contributions for this analysis. 

So from (|22j) , we can write down the branching ratios: 



Br(^ Q -> £p£plp) = lf a Br(4 tpP a Vp) , 

1 {£ a -> £pv a vp) 

~ |A (3Q | 2 (12sin 8 6» t0 - 12 sin^ 9 W + 3 sin 4 6 W ) Br(^ Q -> £pv a pp) , for a = fi,r, (40) 

and for a = r only 

Br(£ a -> £Jplp) ~ \X aa \ i {%wcL % e w -8sin 6 6» u) + 2sin 4 6l 1 „)Br(£ Q -> £pv a V fi ) , (41) 

Br(^ Q -> £ a £p£p) ~ 2|A /3CT | 2 |A 0Q | 2 (sin 2 ^ - l/2) 4 Br(^ -> £pv a Vp) , (42) 

where we have kept onl y th e leading order terms. 

Using the data from [241 ], we can derive the following limits for the various elements of A. In (|40| . there are three 
kinematically allowed processes (fj, — > 3e, r — > 3e, r — ¥ 3fi) and one gets 

|A e „| <4.7x 10- 6 , (43) 

I A eT I < 2.1 x 10" 3 , (44) 

(AjutI < 2.0 x 10~ 3 , (45) 

while (|4"Tj) has two possibilities (r — > efj,JI and r — > \iee), yielding 

|A eT | < 3.4 x 10~ 3 , (46) 

lA^I < 2.2 x 10- 3 . (47) 



Finally, we have 



from another two possibilities (t — > e/x/i and r — > /lee) allowed by (|42[) . 

As expected, these LFV processes provide a stronger set of constraints than those derived in ([37| to (|39| from the 
previous section. 



|A Me ||A MT | < 3.6 x 10~ 3 , (48) 
|A e JA er | < 3.3 x 10" 3 , (49) 



5 Note that the LFV branching ratios quoted in [24 ) is in fact the experimental values for Br(Z — > tain) + Br(Z — > C a £p)- Therefore, the 
expression in (136 I t must be multiplied by a factor of 2 before applying the experimental numbers. 
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FIG. 1: Lowest-order diagrams that are relevant for the amplitude calculations of LFV decays (l a — > Ip^f) and anomalous 
magnetic moments of SM leptons (when a = /3) in the unitary gauge. Subscript j denotes the flavor of the internal leptons, 
and is summed over in the computation, (a) The case mediated by Vj corresponds to the usual diagram studied in standard 
electroweak theory, while the L~~ diagram comes from the new interactions; (b) & (c) are new contributing diagrams involving 
the Z boson and physical Higgs, H , respectively. 



V. RADIATIVE i -> £' 7 DECAYS VIA ONE LOOP 

It is clear that given the continual experimental effort on improving the bounds associated with LFV radiative decays 
of charged leptons (£ — > £'7), 6 any new contribtutions to these interactions originating from the exotic doublets must 
not be overlooked. Therefore, in this section, we calculate the effects due to having additional doublet particles L~ 
and L running inside the one-loop diagrams as depicted in Fig. [TJ 

To set our notation, consider the following generic transition amplitude for £ a — > £p"f. 

T(£ a -> £37) = up (A + Bj 5 )ia p „q"e p u a , a pv = i [7^,7,,] /2 , (50) 

where A and B correspond to the transition magnetic and electric dipole form factors 7 , while q v and e p denote the 
photon 4-momentum and polarization respectively. 

Applying the modified coupling matrices given in (|19p to (|21[) . we can easily relate the new physics parameter \ a p 
to these LFV processes. When explicitly computing the diagrams in Fig. [1] in the unitary gauge, 8 our strategy is 
to perform the calculations in terms of the generalized renormalizable gauge H3, and subsequently, taking the 
limit £ — > 00 to obtain the desired results. 9 Moreover, we will work exclusively in the mi j <C Mw,z,h and mp <C 1 
limits (where rug, and mp represent the masses of the internal j-flavor and the final state SM lepton respectively), 
dropping any sub-leading order terms in the process. In these limits, one also finds that amplitudes A and B become 
identical and we may simply pick out the coefficients associated with the up(l ± 7s)(2p • e)u a components in (|50|) . 
where p denotes the momentum of £p, to get our final expressions. 10 

After the dust has settled and to leading order in A, we obtain the following expressions for the amplitudes of the 
various one-loop contributions (superscripts and subscripts denote the type of internal leptons and bosons involved 
respectively): 



A Z-- = 5w^^ rr .) (jtf-lyt) ^ [/,(„,) + Mu,j)] . w, s MllMl , (52) 
A l z = / 9 ° 1 - - sin 2 6 W + - sin 4 W Xp a , (53) 



8^ 2 V2 V 3 "'3 
iGFiri a e \ ^ v 2 



6 See for example the review in [3(J. 

7 It is understood that A and B are dimensionful quantities when written in this form. Also, we have absorbed the extra factor of i into 
B which is usually factored out in the definition of the electric dipole moment term. 

8 Note that these are the only graphs wc need to consider in this gauge. 

9 We have adopted the definition of £ as used in modern textbooks l2g| . [29ll , which is equivalent to the parameter 1/f as appeared in |27|| . 
A more detailed discussion on this procedure can be found in Q|, as well as many standard textbooks, see e.g. [28j ■ 
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A l- _ iG F m 

A H — 



Ai _ iG F m a e ( m 2 

4_ wr°te 



Ml 



8ir'- 



with 



h(x) 
Mx) 

hA{x) 
Hb{x) 
Uix) 

hi?) 

fe(x) 



-10 + 43x - 78x 2 + 49a; 3 - 4x 4 - 18x 3 Inx 

12(x- l) 4 
4 - 15x + 12a; 2 - - 6a; 2 lnx 

2(x - l) 3 ' 
4 - 9a; + 5a; 3 - 6x(2x - 1) In x 

12(a;- l) 4 ' 
2(1 - x 2 + 2xlnx) 

{x- l) 3 ' 
-11a; + 18a; 2 - 9a; 3 + 2a; 4 - 6a; lnx 

48(i- l) 4 ' 
— 3x + Ax 2 — x 3 — 2a; In x 

8(1- l) 3 ' 
— 2x — 3x 2 + 6a; 3 — x 4 — 6x 2 In x 
24(x- l) 4 ' 



(55) 
(56) 

(57) 
(58) 
(59) 
(60) 
(61) 
(62) 
(63) 



In the above, m v . , m Q and denote, respectively, the masses of the j-navor neutrino, the decaying SM lcpton £ a 

and the exotic L particle. Note that (|S"Tj) is nothing but the contribution due to the SM electroweak interactions when 
neutrinos carry a nonzero mass. It is well-known that the size of this is negligible [28|, [3(1 [3l[ as it receives a m 2 ./M' 2 v 
suppression. 

With these amplitudes and the general formula for the total decay rate, 



-in 



Ayy + A^y 



A% + A L Z 



A l H + A L H 



(64) 



we obtain the branching ratio (after dropping A^,A H ~ 0) 



Br(* a -> = ^ 

Z7T 



1 - | sin 2 6 W + t sin 4 6 W ) X Pa + £ £ (?) ^ Mr 2 (yt)^ 



h{wj) + f 2 (wj) 



+f 3 A(z 3 ) + fsB(Zj) + U(Zj) + 3f 5 ( Zj ) + 2f 5 (h 3 ) + f a (hj) 



Br{£ a ^z/ a F^) , (65) 



where a e is the fine-structure constant. Taking Mj ~ (100) GeV (this is around the global lower bound for h eavy 
charged leptons from current experiments [lij]) 11 for all j, and assuming the Higgs mass, Mh, is about 114 GeV [32| . 



Note that the resulting bounds will become less stringent as we increase the value for Mj. Also, we would like to remind the reader 
that the reason we need to specify a size for Mj here is solely for the evaluation of the loop functions d 5 71) to J 63 1) . We have checked 
that the value for these loop functions would only change by a small amount even when we take the very large Mj limit. 
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the experimental limits 24J on Br(/i — > e7), Br(r — ^ e7) and Br(r — > ^7) then lead to 12 

I A ejU I < 2.4 x 10 -5 , (66) 

|A eT | < 3.0 x 10 3 , (67) 

|A Mr | < 3.5 x 10~ 3 . (68) 



Although these bounds are weaker than those displayed in (|43|) to (|45|) . the derived expressions above will be very 
useful when the expected improvement in the experimental bounds are realized in the near future (besides the upper 
limits, (|67[) and (|68| . are only marginally bigger than their counterparts). Currently, the MEG experiment [33| located 
at Paul Scherrer Institute (PSI) is planning to reach a sensitivity of O (l0 -13 ) for the /i — > e^j branching ratio, which 
is a significant improvement compare to the current limit of Br(/i — > ey) < 1.2 x 10 -11 (34]. In addition, the Super 
KEKB project [35| will provide an excellent platform for investigating LFV r decays at an unprecedented precision. 
As a result, the bounds on r — > and r — > ^7 are also expected to tighten, providing a stronger constraint for all 
off-diagonal couplings X a p than presented here. 



VI. fi-e CONVERSION IN ATOMIC NUCLEI 



Muon-to-electron conversion in muonic atoms provides another excellent testing field for tree- level FCNC. This is 
because coherent contribution of all nuclcons in the nucleus can enhance the experimental signals, and hence the 
fX-e-Z vertex may be probed at great precision. Given that this is the same vertex as appeared in the loop graphs in 
Fig [lb, the test for pi-e conversion, therefore, plays a complementary role to the investigation of fi — > e~f in the probe 
for physics beyond the SM as the two processes are induced differently. 

In what follows, we shall assume that the only contribution to the fi-e conversion rate in our setup comes from 
exchanges with the Z bosons. This approximation is sensible because the cases mediated by the photon and the Higgs 
are suppressed by loop effects and M^ 1 respectively. So, assuming only SM interactions operate in the quark sector, 
we obtain the following effective interaction Lagrangian for the /U-e transition (after integrating out Mz): 

= G F l e "f v {kv ~ kAjs^/i [q u Ju( v u + aulb)qu + q,ilv( v d + a<275)<?d] , (69) 

where q Utd denotes the u, <i-quark field while 

/l o \ 1 1 4 2 12, 

k v = -k A = I - - sin 9 W J A eAl , a u = -a d = -- , v u = - - - sin 8 W , v d = - - + - sin 6 W . (70) 

Appealing to the general result obtained from FCNC analysis with massive gauge bosons in (2f|, the branching 
ratio for fi-e conversion in nuclei (for nuclei with less than about 100 nuclcons) is found to be 

where p' e (E' e ) is the momentum (energy) of the electron, r^ p represents the total nuclear muon capture rate for 
element A, and Z (Z e g) is the (effective) atomic number of the element under investigation. In (f7"Tj) . F(q' 2 ) is the 
nuclear form factor which may be measured from electron scattering experiments [36| while 

Q = (2Z + AO v u + {Z + 2N) v d , (72) 

with Af denoting the number of neutrons in the nuclei. 

Given that one of the best upper limit on the fj,-e conversion branching ratio is obtained from measurements with 
titanium-48 (||Ti) in the SINDRUM II experiments [13 : 

„,w„ r(/i~ Ti — > e~ Ti) 1Q 

B^ e = ^ < 4.3 x 10- 12 , (73) 

cap 



If we compare this set of limits with the corresponding ones in 0], we notice that these bounds are somewhat stronger. The reason for 
this comes from the fact that Q omitted the additional W-mediated graph where the SM leptons couple to an internal doubly charged 
exotic triplet. As a result, the small accidental cancellations in the numerics (as happened in these numbers here) did not happen there. 
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we shall use the parameters for element ^Ti in (fTTj) to deduce our bound. 13 Following the approximation as applied 
in [H, we take p' e ~ E' e ~ m*, and F(q' 2 ~ -m£) ~ 0.54. In addition, we have Z off ~ 17.6 for ||Ti [H and 
r^ p = r^ p ~ 2.59 x 10 6 s" 1 gl| . Hence, ([711) and (£3) combine to give 

|A ep | < 1.9 x 10" 6 . (74) 

As hinted earlier, the bound displayed in (|74[) is indeed the most stringent one on \X efJ ,\- Moreover, given that new 
/i-e experiments are being planned, respectively, at J-PARC and Fermilab by the COMET (and PRISM/PRIME) [HJ 
and Mu2e [43| collaborations, this bound is expected to be further strengthened in the near future. 



VII. GLOBAL FIT ON THE ELEMENTS OF A AND SOME CONSEQUENCES 

In this section, we bring together all the results obtained thus far and perform a global analysis on the elements of 
the A matrix, which are key to determining the new physics effects from the doublet leptons, L. For convenience, a 
summary of all constraints derived in the previous four sections are collected in Table [II] 



parameter(s) 


process 


limit on BR 


constraint on A's 


A ee 


Z -¥ e"e + 


3.363 ±0.004 x 10 -2 




'W 


Z — > /-L /-L + 


3.366 ±0.007 x 10~ 2 


< io~ 3 


A T T 


Z -> T--T + 


3.369 ±0.008 x 10" 2 






Z e ± f^ 


< 1.7 x 10~ 6 


< 6.6 x 10~ 3 


|Ae M | 


fj,~ — > e~e~e + 


< 1.0 x 10~ 12 


< 4.7 x 10~ 6 


/i — > cy 


< 1.2 x 10 -11 


< 2.4 x 10" 5 




H~e conversion 


< 4.3 x 10" 12 (Ti) 


< 1.9 x 10" 6 




Z e ± r =F 


< 9.8 x 10~ 6 


< 1.6 x 10 -2 


| Aer | 


r~ — > e~e~e + 


< 3.6 x 10" 8 


< 2.1 x 10" 3 


t~ — > e~/i~/i + 


< 3.7 x 10~ 8 


< 3.4 x 10~ 3 




r — 5- ej 


< 3.3 x 10~ 8 


< 3.0 x 10" 3 




Z -> 


< 1.2 x 10~ 5 


< 1.8 x 10" 2 


|A, 1T 


r~ — >• fi~ 


< 3.2 x 10~ 8 


< 2.0 x 10~ 3 


t~ — > fi~ e~ e + 


< 2.7 x 10 -8 


< 2.2 x 10 -3 




T — > /i7 


< 4.4 x 10" 8 


< 3.5 x 10~ 3 


| A^ e | | A^r | 


r — i> e + fi~ (i~ 


< 2.3 x 10" 8 


< 3.6 x 10" 3 


| A eAi | A er | 


t — > fi + e~e~ 


< 2.0 x 10 -8 


< 3.3 x 10" 3 



TABLE II: A collection of all constraints on the elements of A = v 2 YM~ 2 Y i ' /2 from processes studied in the previous four 
sections. 



Studying the results listed in Table [TT1 and recalling that \\ a p\ 
overall fit for the elements of A: 



|A^q,|, it is not difficult to obtain the following 



(\^ee\ 


|A e/1 | 


|A 


|A M e| 


l-Wl 


|A 




|Ar^| 


|A 



< 



10~ 3 1.9 x 10~ 6 2.1 x 10~ 3 
1.9 x 1(T 6 10~ 3 2.0xl0~ 3 
\2.1 x 10~ 3 2.0 x 10" 3 10" 3 

Result ([75)) is one of our major results in this work. Using the definition, A = v 2 Y M~ 2 Y^ /2, and taking M 
O (100) GeV (for all flavors), a rough estimate of the size for the new couplings may be obtained: 



(75) 



Yi 



<O(l0" 2 ) toO(lO- 3 ) , for all 



i,3 



(76) 



Although the value quoted in the experiments with gold (Au): r(/i — Au — > e _ Au)/r*° < 7 X 10 — 13 f38| is smaller than the one in 
1 1731 . theoretical calculations |39( have shown that for very heavy elements (atomic number Z > 60) like Au, the fi-e conversion rate 
is actually suppressed. Therefore, this does not necessarily indicate a better bound on the rate, especially when the estimation of the 
nuclear matrix element for such heavy nuclei can carry large uncertainties. 
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If the exotic particle mass M is heavier than the value used above, the upper limit for \Yij\ will be increased. 

Furthermore, if we assume that the new flavor changing physics due to the presence of these exotic L's are the only 
source of LFV, one may derive model-specific bounds on the various processes predicted by this model from studying 
the ratio between the different branching ratios: 

Br(Z — » e ± ^i T ) ~ 4.8 x 10 _1 Br(/i-e conversion in Ti) , 
Br(/i~ — > e"e~e + ) ~ 3.9 x 10~ 2 Br(/i-e conversion in Ti) , 

Br(/z -> cy) ~ 1.9 x 1CT 2 Br(/i-e conversion in Ti) , (77) 

for the processes involving |A eM | 2 . Whereas for |A eT | 2 and |A Mr | 2 , one may write 



and 



Bt(Z -> e ± r T ) 


~ 7.0 x 


10 1 


Br(r~ - 


e~e"e + ) , 




Br(r _ — > e~//~/i + ) 


~ 6.5 x 


10" 


1 Br(r" 


— > e~e"e + ) , 




Br(r — > ey) 


~ 4.6 x 


10" 


1 Br(r" 


— > e~e~ e + ) , 


(78) 


Bi{Z -> ^r^) 


~ 7.1 x 


10 1 


Br(r" - 


-> /Oi"m + ) ' 




Br(T~ — » /i"e~e + ) 


~ 6.9 x 


10" 


1 Br(r" 






Br(r — > /i7) 


~ 4.6 x 


10" 


1 Br(r" 




(79) 



respectively. Whenever required in the above, we have used M ~ O (100) GeV [13 - 14 The main point is that we can 
translate these into model-specific bounds for certain LFV processes which may be used to falsify this theory. For 
instance, applying the experimental limits on the right-hand side of (|77j) implies that this model demands Br(fi — > 
e l) ^ 8.2 x 10" 14 . Observe that this is a couple of orders stronger than the limit set by current experiments. As a 
result, a future detection of this LFV process above this rate will invalidate the predictions of this minimal extension 
to the SM, and point to the existence of other new physics in the lepton sector. Similar conclusions may also be 
drawn from other processes displayed above. 



VIII. CONTRIBUTION TO LEPTON ANOMALOUS MAGNETIC MOMENT 

While in Dirac theory the gyromagnetic ratio of a spin-1/2 particle is predicted to have a value of g dlrac = 2, it 
is well-known that quantum field theory gives a correction to this number via loop effects. The deviation from the 
Dirac result of 2 is usually parameterized by the dimensionless quantity (a denotes the flavor) 

a a = — , where g a is the actual value of the gyromagnetic ratio, (80) 

known as the anomalous magnetic moment. It is related to the lepton magnetic dipole moment \± a = —e(l+a a )/2m a s, 
where s is the unit spin vector. In terms of the parameters from quantum field theory, a a = -F^T = 0), when the 
form factor expansion for a general lepton-photon amplitude is written as 

T{t a ^l' crf ) = -ieu' c 

where q" is again the photon momentum (see (|50[) for notations). 15 Therefore, the precise contribution to a a from the 
SM (and indeed any other theories) can be calculated by considering all the relevant loop diagrams for the i 7 2(0)-term. 

While the anomalous magnetic moment for the electron, muon and tauon can all be very important in their own 
rights, given the present experimental and theoretical development, a M is the most interesting observable to examine. 
This is because when combining the fact that significant contributions to the overall predicted a M value come from 



2m n 



2m„ 



■ CV75 q 



e H u n 



e > 



(81) 



We have checked that taking a larger value for M would only change the numerical results by a small amount. 
Note that the lepton electric dipole moment is proportional to F^q 2 = 0). 
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every major sector (QED, electroweak, hadronic) of the SM j4J:Eal with the ability to experimentally measure to 
extremely high accuracy [46|, |4?J , the SM as a whole can be scrutinized, and any discrepancies between theory and 
experiment would be a strong indication of new physics. On the other hand, although a e have been measured to 
extraordinary precision (hence providing a very stringent test on QED and the value of the fine-structure constant a e 
13 IHI): hs low sensitivity to the contributions from strong and electroweak processes means that any hypothetical 
modifications to these sectors (due to new physics) would not be easily detectable. As far as a T is concerned, even 
though its much heavier mass would in theory imply better sensitivity to any new physics than a M , its usefulness 
has been limited by the relatively poor experimental bounds. In fact, the best current limits set by the DELPHI 
experiments [5(| are still too coarse to even check the first significant figure of a T from theoretical calculations. 
Currently, the experimental values for a e [48|, [13] and a T [5(| are given by 

af xp = 115 965 218 073(28) x 10" 14 , (82) 
a^ p = 116 592 089(63) x 10~ n , (83) 

«?* =[ <1 - 3X ^ a (84) 
| > -5.2 x 10- 2 . V ; 

Focusing on the muon case, one finds that the discrepancy between experiment and the SM estimate is about 4.0a 

Aa M = a^ x P - af* = 316(79) x 10" 11 . (85) 

If this difference is real (rather than caused by incorrect leading-order hadronic approximation 17 ), then there must 
be some new physics at play. In the following, we investigate whether the presence of the exotic doublets can affect 
this quantity in a significant way. 

Calculating the anomalous magnetic moment using the modified electroweak couplings of (fTTj)) to (|2"Tj) is in fact 
analogous to the computation for LFV l a — > ipj done in Sec. [V] The only differences are that a = j3 here and we do 
not take the zero mass limit for the final state lepton. Otherwise, the three main types of one-loop diagrams we need 
to consider arc as depicted in Fig. [TJ Working in the unitary gauge again, and noting that for the magnetic moment, 
the part associated with 75 in the general amplitude 

T(t a -> Cl) = < (C + £> 75 ) iG pv q v e ( 'u a , (86) 

is not needed. Hence, in the computation, we pick out the terms that are proportional to u a (2p- e)u a , where p is again 
the momentum of the incoming t a . Employing a similar notation system as before, the amplitudes of the one- loop 
diagrams from Fig. [T]for the case a = /3 are given by (to leading order) 

iG F m a e f 5\ 

c t- = ! gm^ £ £ (YM~') aj (M- 1 ^)^ [f 7 ( Wj ) + 3/ 8 K) + h{w 3 ) - 1] , «y ee M^/M^ , 
iGFm a e (2^ : n _. 2 4 2\q^q, ^ - 4 /1 



z = is^jt U (1 + 2 9w ~ 4 ° w) + ° w ~ 8 9w + 3) J ' (89) 

cf = y (™ l ) aj (^~ 1?t ) JQ [2/3A(Zj) + 2/3sfe) + /4( ^ } + 2/sfe)] ' z > = ^ /M| ' (90) 



C e H = (1 - 4A QQ ) O {ml/Ml) ~ , (91) 

°" = iS^K X> Y Y a] Mj>Yl [4f 5 ( hj ) + 2f 6 ( hj )] , hj = Mj/Mfj , (92) 



Recently, this discrepancy was re-evaluated by the group in |51| and found to be about 3.3<r only: Aa M = 261(80) X 10 — 11 . 

Although this possibility is not completely ruled out, shifting the hadronic cross-section to bridge this gap will naturally increase the 

tension with the lower bound on the Higgs mass, both from LEP \Sj and the SM vacuum stability requirement |4T| . 
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where /3a(x) to fe(x) are given in (f59|) to (|63|) and 



7-33a; + 57a; 2 -31a: 3 + 63; 2 (3a;-l)lna! 
M*) = ^31)4 > (93) 



1 - 4a; + 3x 2 - 2a; 2 In x 



M*) = T— TT3 • ( 94 ) 



, 3 - 10a; + 21a; 2 - 18a; 3 + 4a; 4 + 6a; 2 In x , s 

M*) = ^3Tji ■ (95) 

Comparing (|86[) with the form factor expansion of (|81|) . the anomalous magnetic moment can be written in terms of 
the amplitudes computed above: 

a a = F 2 (0) = (C^r + C&~~ + C| + Cf + C H + Cl~) . (96) 

Note that the result given in (l96l) contains the usual SM electroweak component of a a , as well as the contribution 
induced by the new physics. Examining our results, we see that the terms which are not proportional to X aa in 



([57)) and (|89p sum up to give the usual prediction of the anomalous magnetic moment from the SM [27J. Removing 
this component from (|96[) and using the values for X aa given in Table UH we obtain the following estimate for the 
anomalous magnetic moments coming from the new physics associated with the exotic L particles: 

Aa^ ~ 6.6 x 1(T 16 , (97) 
Aa£ ~ 3.2 x 10~ n , (98) 
Aa^ ~ 1.1 x 10" 8 , (99) 

where we have again assumed Mj ~ O (100) GeV and Mb ~ 114 GeV. Looking at the results (|i?Tj) to (|9"9")l , we see 
that these contributions are at least one order of magnitude less than the experimental errors given for the quantities 
listed in ([82]) to (|84|) . Therefore, these exotic L's cannot help to explain the muon g — 2 anomaly nor can their effects 
be easily distinguishable from the SM components in these experiments. 



IX. COMMENTS ON COLLIDER SIGNATURES 

If the exotic doublets L do exist, it would give rise to collider signals which may be detectable at the LHC. Given 
that the exotic mass M is not too massive and assuming some favorable conditions are met, such signatures can be 
quite distinctive as pointed out recently in Q . In this section, we summarize some of the key features one can expect 
from the presence of these exotic particles. 

To make the situation more clear-cut, let suppose Yukawa couplings Y are small enough such that L production 
is predominantly mediated by SM gauge interactions. The relevant production mechanisms at the partonic level are 
then given by: 

qq^j,Z^L-L+ ,L—L++ ; ud -> W + -> Z"Z++ . (100) 

With the planned energy of y/s = 14 TeV and luminosity at 300 fb~ 1 yr _1 at the LHC, it is estimated that the pair 
production cross-section from pp collisions ranges from a(pp — > LL) ~ a few fb for M ~ 500 GeV to just under 10 3 fb 
for M ~ 200 GeV [|[ (assuming no cuts arc made). Once produced, the singly and doubly charged exotic leptons 
may decay into SM particles via the interaction terms depicted in (|T6|) to (| L8|) : 

L-^rZ,rH and L~ -> i~W~ . (101) 

Note that the decay mode L~ — > vW~ is suppressed. Furthermore, we shall assume that the Y mediated processes 
dominate over the electroweak decay L — > L~tt~ , which is, strictly speaking, allowed because of the mass splitting 
of the doublet components due to quantum corrections (52| . As a result, we may treat both L~ and L on an equal 
footing. Consequently, we see that the three relevant states fcovapp collisions (L~L + , L~~L ++ and L~L ++ ) will lead 
to a generic £+£-XX' signal, where XX' denotes one of ZZ, HH, ZH, W + W~,ZW+ or HW+ . 
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For instance, the casepp — s- 7, Z —> L~~ L ++ will give rise to a £ + W + £~W~ state. Subsequently, the W bosons may 
decay leptonically or hadronically leading to one of the final states (in order of descending branching fraction): £ + £~j, 
tr 1 l ± $j'£f'j or £ + £~£ + £~$ T , where an d j denote the generic missing transverse energy and jets respectively The 
relevant SM background processes for this case are pp — > iiW , ZW + W~ (and perhaps pp — > it). Amongst the three 
final states, £ ± £ ± $ T £ T j seems to be the most promising for probing the mass of M as one may gain information 
from studying either the two same-sign leptons invariant mass distribution or the invariant mass of the opposite-sign 
lepton with two jets. 

Another interesting case to consider is the pp — > W + — >• L~L ++ decay chain. Although there are more possibilities 
for the intermediate state £ + £~XX', if one assumes that any subsequent W boson decays only leptonically, then the 
final signal one gets is either £ + £ + £~ $ T or £ + £ + £~ $ T j . The typical SM background one must confront with here is 
coming from vector boson decays (e.g. from pp — > ZW + ). According to a recent analysis in |3|, such channels can 
provide another promising way to search for these exotic particles. Note, however, that the cleaner state £ + £ + £~ \$ T 
is only possible if M is low enough. 

X. CONCLUSION 

Given that the discovery of nonzero neutrino masses demands an extension to the lepton sector, it is natural to ask 
what might be the simplest ways that new physics can couple to the SM particles. Concentrating on the lepton sector 
only, we have followed the approach of Q and introduced exotic particles into the SM via some "minimal couplings" . 

Since many of the exotic particles defined in this way turn out to be equivalent to those studied in other new 
physics models, we identify that only a couple of possibilities remain unexplored in the literature. While the work of 
[2| focused on one of them, in this paper, we have presented the analysis for the last remaining possibility, namely, 
the exotic doublets L which couples to the RH charged lepton singlet. 

Using a formalism similar to that presented in we have defined the key quantity, A, which encapsulates the new 
physics effects caused by the introduction of these exotic L's. In particular, we note that the off-diagonal entries of 
this A matrix are the origins of any new FCNC phenomenologies. By invoking the limits from low-energy experiments, 
constraints are then placed on these entries that control the coupling strength to the exotic L's. Such an investigation 
can be quite useful given that these minimally coupled particles may give rise to definite collider signatures at the 
LHC in the future [3|. 

In this paper, the processes considered include leptonic Z decays, LFV £ —> M, £ — > £'7 decays, as well as fi-e 
conversion in titanium nuclei. The collection of constraints are displayed in Table HU We have found that diagonal 
elements, |A aa |, could be as big as O (l0~ 3 ) only, while the most stringent bounds for the off-diagonal values are 
from LFV £ — > 3£ decays except for |A eM | which has its strongest bound coming from the p-e conversion process. We 
anticipate that some of these limits will improve significantly when the next generation of experiments have reached 
their proposed sensitivity. 

Finally, the contribution to the lepton anomalous magnetic moment in this model is calculated. The explicit 
computation of the relevant lowest-order loop graphs shows that any potential contributions is far too small to be 
detected in experiments at the present time. As a result, introducing this type of exotic doublet of leptons into the 
SM cannot resolve the muon g — 2 anomaly. 
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